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\ Abstract 

We analyze the model of topological fermions, where charged fermions 
qq ' are treated as topological solitons. We discuss vibrations of soliton 

■ shapes. It is shown that depending on the power of the potential term 

(discrete parameter m) of the model Lagrangian the spectrum of normal 
mode frequencies can be discrete (for m = 1) or continuous (for m > 2). 
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1 Introduction 

The well-known success of the Skyrme model to the description of short-range 
forces and properties of strongly coupled particles makes it worthwhile to ex- 
tend the "skyrmion philosophy" to long-range forces and physics of electrically 
coupled particles. The so-called model of topological fermions (MTF) [21 [3] 
proposes a realization of such an idea. 

The model has three independent degrees of freedom parameterizing an 
SU(2) field 

Q(x) = coso(x) + ian(x) sin a(x), (1) 
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where o are the Pauli matrices and n(x) is a three-dimensional unit vector in 
internal ("colour") space 0. Due to the constraint = 1 this vector has 

two independent degrees of freedom. Both fields, a(x) and n(x), are functions 
of the Minkowski coordinates x^ = (ct, x, y, z). 
The Lagrangian density of the MTF reads 

£ = -^(^-^ + A(go)), (2) 

where W v is the curvature tensor 

BT = f> A f u , with the connection f M = —Tr( ad^QQ^). (3) 

2i 

The potential term is given by 

1 /TrO\ 2m 1 
A(9o) = ihr =-cos 2m «(x), m = 1,2,3,... (4) 
r o V A J r o 

The model contains two parameters, the fine-structure constant, a/, and a 
dimensional parameter ro- 

Note that the "curvature term" —jR^ ■ R^ u is proportional to the Skyrme 
term, but the so-called kinetic term of the Skyrme model does not enter the 
Lagrangian density (T2J) in order to allow for electromagnetic fields and forces 

Due to its Lagrangian density the MTF has different properties than the 
Skyrme model at r — > oo [2], |3] . In the Skyrme model the chiral field U ap- 
proaches the trivial configuration, U — > 1. In the MTF the field configuration 
for r — > oo is determined by the minima of the potential characterized by 
ac(x) = | and arbitrary direction of n, 

Q(x) = ian(x) at r — ► oo. (5) 

As a result the Q field can form a hedgehog configuration and the field a(x) 
describes the profile of a charged soliton with properties of an electron, whereas 
the field n(x) is related to the dual electromagnetic field strength [21 [3] by 

*fnu{x) = ~ , 6 ° [d^n(x) A d v n{x)] ■ n(x). (6) 



^^We use the summation convention that any capital latin index that is repeated in a 
product is automatically summed from 1 to 3. The arrows on variables in the internal 
"colour" space indicate the set of 3 elements q = (91,92,93) or = (01, 02, 03) and qa = 
Qk^k- We use the wedge symbol A for the external product between colour vectors (qAa)A = 
zabcIbvc- F° r the components of vectors in physical space x = (x,y, z) we employ small 
latin indices, k and a summation convention over doubled indices, e.g. (E x B)j = 
CijkEjBk- Further we use the metric 77 = diag(l, — 1, — 1, —1) in Minkowski space. 
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The field strength reads = --e^* f p ° with e mM = 1. 

The model has two types of excitations, which are related to n and a 
degrees of freedom and have different physical meaning. 

One type (connected with fluctuations of the field n(x)) is realized as elec- 
tromagnetic field evolved by the charge in the wave zone (where a — > |). 
Some properties of such fluctuations were already studied, appropriate classi- 
cal equations of motion for the field n were derived [3] and explicit solutions 
of these equations of motion, which behave like electromagnetic waves, were 
found pp. 

The subject of the present paper is to study another type of MTF exci- 
tations which is generated by fluctuations of the field a(x) and realized as 
vibrations of the soliton shape. 

The paper is organized as follows. In Section [2] we derive the Lagrangian 
density for small a-fluctuations. Than, in Section[3], we calculate normal modes 
of the fluctuations. We find that the spectrum of the mode frequencies is very 
different for m = 1 and m ^ 2, where m is the power in the potential term 
(1211) . The conclusions are given in Section HI 



2 a-fluctuation in the MTF Lagrangian den- 
sity 

We will start from the second order variation terms of the MTF Lagrangian 
density 



5C = -^5(\R^W \ 

a the 



4-7T 



C 2 q d qo A - UCQ) 2 d 2 q A - (d.C A d u ( )BT + (d,C A fl (r„ A RT) 

+1(8^ a &)(d„c a n - a f,)(^c a n , 

where ( = C(^ x ) are three independent variation parameters introduced by 
Q^Q' = e^Q = 

i ,_n\ . _x T* / . __>_.\ 

(8) 



(1-|C 2 )+^C (qo+iqv 



+ io 



;i-K 2 )?+?oC + ?ac 



— * 

In Eq. (J7J) terms linear in ( vanish due to the equation of motion. 

To simplify calculations we introduce in (J7J), (jSJ) and later on the notations 

qo{x) = cosa(x) and q(x) = n(x) sina(a;) and use spherical coordinates 9 and 

in colour space. 
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We will also use the spherical coordinates, r, (p for the physical space and 
specify the hedgehog soliton by n = -. One gets 



n. 



d 9 n 
sin?? 



By rotation of eg and eg and e^ with angle a we get 

= cos a eg — sin a e^, e v = sin a eg + cos a i 

with 



(eg A e v ) = n 



and 



(c^n A c^n) • n = sm-d. 



(9) 



(10) 



In spherical coordinates the covariant and contravariant components of the 
connection and the curvature tensor (ED read 



R, 



a [r) n, 

sin a eg, 
: sin •d sin a e v , 
= sin t? sin 2 a ft, 
= sin^sinaa'eg, 



-a r n, 



sin a a' e 



p? _ _ l 
fV = _ 

R T ^ 



-s sin a eg 
i 



- 2 ■ a sin ae„ 

r z sm # '/ 
sin a 
r 4 sin i9 ' 
sin a a' 



r sm 
sin a a 



(12) 



^eg, 

2 ^rj- 



After some algebra one arrives at 



5C 



ajhc 
An 



2 sin a a 



r 2 sin # 



|C 2 cos a d qo A - \C sin 2 a d 2 qo A - f^f A d v ( 
'(^CA9 r C)-eg 



n 



2 sin a a' 



(9 r CA^C) -e; + 
+ ^(«' 2 + ^)P,C A C) eg + -^(«Vf A C) e„] 
+|{a'(9 r C A f») + ^[(^C A eg) + -^(^C A e^)]} 5 



2 sin 2 a a' 



{drC A C) n 



i^ r» / ^ — *\ / r\n », — #\ 1 / /D e in 2 

r 



sin" o / 



(<9 M C ■ n)(^C ■ fO + |(«' 2 + ^)[(^C eg)(^C eg) + (<9 M C 3,)(0"C 3,)] 

(13) 



Now let us consider "a-fluctuations" , which correspond to the following choice 
of the parameter ( 

C(t,x) = <P{t,r)n, ft =^7. (14) 



Then the fluctuation of the Lagrangian density reads 

oifhc 



5C 



An 



+ 



sin 2ct 



□a + — cos 2m a ( (2m - 1) tan 2 a - 1 



(15) 
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Using the equation of motion 



m 



cos 2m a 



sin a cos a sina. _„ cos a _ _„ sin a 
2 + d^a — d^ad^a — Da 



(16) 



the expression (Tl5l) reduces to 



5C 



ctfhc fsin a 



Aix 



d, 



sin a(l + 3cos2a;) sin 2a 



2r 4 



<9„r<9 M a 



sin 2 a 



_ sin 2a ^ m 

i) l ,<\()"<\ + 



, Da + ^cos a(2m-l)tan 2 a 
2r l 



(17) 



Introducing a new field variable 



a/2 sin a , 
$ = cf) 



one arrives at the final expression for the fluctuating Lagrangian density 
4tt {2 3>m 



5C= °^{ld,M^-& 



l + 3cos2a m(2m — l)r 2 9m 9 
■ H „ „ cos 2m ~ 2 ct 



2r 4 



(19) 



(20) 



3 Normal modes of fluctuations 
3.1 Fluctuation potential 

From (fl9~I) we get the Lagrangian for the a-fluctuations. It has the standard 
form 

L[$] =T[$] -C/[$], T[$] = |fd 3 x<i> 2 
[/[$] = y d 3 a; [|(V$) 2 + $ 2 V(r1 

= y d 3 x $ [-±A + V(r)] $, 

_ 1 + 3cos2a m(2m - l)r 2 2m _ 2 _ tKp) 
1 ' ~ 4r 2 + 2r 4 C ° S " ~ r 2 ' 

where p = r/r$. 

To calculate the potential V(r) one needs to know the profile function 
a(r). The solutions for the profile function with different m are discussed in 
Appendix [A] The behaviour of cosa(r) for m =1,2,3 and 4 is displayed in 
Figure [Q 

From Eq. ( 1201) one immediately learns that for m = 1 the fluctuation 
potential increases as V(r) ~ r 2 at r — > oo, while for m > 2 it decreases as 
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Figure 1: At the left panel: cosa(r) for m = 1 (solid line), m = 2 (long-dashed 
line), m = 3 (short-dashed line) and m = 4 (dot-dashed line). At the right 
panel: the potential v in dependence on p = r/r . 



V(r) ~ r -2 , see Figured] Outside the region 0.3 < p < 2 the potential is very 
similar for all m > 2. 

For m =2 and 3 we can give explicit expressions for the potential 



v(p) 



2 2 - 2p z + lOp 
7 ' p 2 (l+p 2 ) 2 

2 + p 2 + 14p 4 
I p 2 (l + p 2 ) 2 ' 



m 



m 



(21) 



where p is defined in 



3.2 Normal modes 

To find the normal modes of the fluctuations we expand r) in an ortho nor- 
mal and complete set of functions rji(r) (see, e.g, [1]) 



$(t,r) = ^ Q (t)^(r). 



(22) 



These functions r]i(r) are determined by the Schrodinger-type equation 

[-iA + V(r)]r 1i (v) = n 2 i r 1i (v). (23) 

Because the potential is very different for m = 1 and m > 2 let us consider 
these two cases separately. 

• m = 1 



There is an infinite number of bound states. Separating angular and 
radial coordinates 77* (r) = ^R n i(r)Yu 3 (r/r) in Eq. (1231) we get for the 
radial wave function 



1 

2' 



~ ~ R nl + 



1(1 + 1] 

2r 2 



+ V(r) 



Rnl — ^nl-^nh 



(24) 
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where 



and a(r) is given by ([35 



1 + 3 cos2a(r) r 2 v(p) 



(25) 



From (|24p and (I25p it follows that for I ^> 1 the spectrum reduces to the 
spectrum of the three dimensional harmonic oscillator. 

At r — > Eq. (ED becomes 



- 2 R ni + 



/(/ + !) + 2 
2r 2 



-Rn/ — 



(26) 



with the solution 



R»i~re, i= l - + J- A + l(l + l) + 2. 



(27) 



At r — * oo it reduces to the asymptotic oscillator equation 

o 1 



- \Ki + \^r 2 R nl = 0, 



K = U = — 



:28) 



with the solution 



2R nl ~ r CT e- 



(29) 



The numerical solutions for the radial eigenfunctions are shown in Fig- 
ure [2] for different angular momentum quantum number I and radial 
quantum number n. So the solutions of the Schrodinger equation ([23]) 
are characterized by three quantum numbers, % = {n,l,l^). The eigen- 
states are obviously degenerate with respect to l 3 . 

In Figure [3] the eigenvalues are compared with the spectrum fl^ of 
the three dimensional oscillator. 

m > 2 

For m > 2 the potential V(r) is repulsive everywhere and there are no 
bound states. At r — > oo the solutions of Eq. (123]) are spherical waves 



»7* 



and 



a;. 



U2 
2 K . 



(30) 



The spectrum is continuous. 
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4 Conclusions 

We derive the Lagrangian for small fluctuations of the soliton profile around 
the hedgehog solution and discuss shape vibrations of the topological fermions. 
It is shown that at r <C r$ the potential term of the Lagrangian is repulsive and 
very similar for all values of the discrete parameter m of the model. At r — ► oo 
its behaviour differs for m — 1 and m > 2. It grows as r 2 for m — 1, while for 
m > 2 it decreases as r~ 2 . As a result the spectrum of normal frequencies is 
discrete for m — 1 and continuous for m > 2. 
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A Profile function 

The profile function is determined from the differential equation 

d 2 q (1 - <?o)<?o 2 2m-i n f*i \ 

■jY + — 2 mp q =0 (31) 

supplemented by the following boundary conditions 

a(0) = 0, a(oo) = |. (32) 

In (13*11) one sets go( r ) = cosa(r) and A(g ) = <?o m ( r )- 
At short distances, p 1, Eq. ( I3TT) is fulfilled by 

g (r) « 1 - /tp 2 (33) 

with arbitrary k. The parameter k is determined by the condition qo(oo) = 0. 

Now let us discuss the behaviour of go at r — >• oo. This behaviour is very 
different for m = 1 and m > 2. Thus we will consider these two cases sep- 
arately. Note, that for m = 2 and 3 there are exact analytical solutions, for 
other m only approximate solutions do exist. 

For m = 1 one gets the asymptotic equation 



d 2 q ., 



dp 2 



p q = at p — > oo (34) 



with solution g ~ p ?e 2 p2 . To connect the solutions (1331) and smoothly 
one can use the following trial function 
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The variation parameter k,q = 0.206796 is determined from the minimum of 
the energy functional 



H[q ] = I dp 
lo 



2\2 



:i - g o 2 ) 

2p 2 



+ (d p q ) 2 + P 2 qT 



For m > 2 Eq. ( !3~Tj) is reduced to the equation 

d 2( l0 . QO 2 2m-l n + 

-— - + — - mp g =0 at p -> oo, 

dp 2 p 2 

which is different from (1340 and has the asymptotic solution 



(36) 



(37) 



g = Ap where A 
For m = 2 



m 2 + 3 
m(m — 1) : 



2(m-l) 2 

and £ 



m — 1 



air 



cos a 



arctan(pi/2 + p 2 ), 

1 pV^+W 



1+p 2 



sm a 



1 + P 



For m = 3 



/ \ / \ 1 P 

a(r) = arctan(p), cos a = — ^ = , sma 



1 + p 2 



aA+? 



(38) 



(39) 



(40) 



see P]. 

For m > 4 one can connect the two solutions ( 1331) and ( 1381) by the following 
trial function 

g = (1 + ^iP 2 + K 2 p 4 )-«/ 4 (41) 



with 



k 2 = ^4 « 



m(m — l) z 
m 2 + 3 



(42) 



and K\ > is a variational parameter. 
For m = 4, 



cosa r 



{/l + Kl p 2 + K 2 p 4 



with ki = 2.98428, k 2 = f§. 
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Figure 2: The radial wave functions R n i for m = 1. The radial wave functions 
are shifted by eigenvalues f^. For comparison we show the sum of the potential 
and the centrifugal energy ■ 
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Figure 3: The spectrum of Eq. ff23l for m — 1 (solid lines). For comparison 
we also show the spectrum Q^i = ^(2^ + / — \) of three-dimensional oscillator 
potential (dashed lines). 



